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Abstract 

The local composite gluon-ghost operator (-j^A^A® + ac a c a ) is analysed in the frame- 
work of the algebraic renormalization in SU(N) Yang-Mills theories in the Landau, Curci- 
Ferrari and maximal abclian gauges. We show, to all orders of perturbation theory, that this 
operator is multiplicatively renormalizable. Furthermore, its anomalous dimension is not an 
independent parameter of the theory, being given by a general expression valid in all these 
gauges. We also verify the relations we obtain for the operator anomalous dimensions by 
explicit 3-loop calculations in the MS scheme for the Curci-Ferrari gauge. 
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1 Introduction 



Vacuum condensates are believed to play an important role in the understanding of the non- 
perturbative dynamics of Yang-Mills theories. In particular, much effort has been devoted to 
the study of condensates of dimension two built up with gluons and ghosts. For instance, the 
relevance of the pure gluon condensate (^A a ^A®S in the Landau gauge has been discussed from 

lattice simulations as well as from a phenomenological point of view |2j . That the operator A 
has a special meaning in the Landau gauge follows by observing that, due to the transversality 
condition d IM A a>M = 0, the integrated operator (VT)~ 1 j d 4 x A a ^A a ^ is gauge invariant, with VT 

denoting the space-time volume. An effective potential for (a^A^S has been constructed in [3], 
showing that the vacuum of Yang-Mills favors a nonvanishing value for this condensate, which 
gives rise to a dynamical mass generation for the gluons. 

The operator A 2 ^ in the Landau gauge can be generalized to other gauges such as the Curci- 
Ferrari gauge and maximal abelian gauge, (MAG). Indeed, as was shown in [I], the mixed 
gluon-ghost operator 1 O = (hA^A® + ac a c a ^j turns out to be BRST invariant on-shell, where 
a is the gauge parameter. Also, the Curci-Ferrari gauge has the Landau gauge, a = 0, as a 
special case. Thus, the gluon-ghost condensate [ \ I^A a ^A a ^j + a (c a c a )^j might be suitable for 
the description of dynamical mass generation in these gauges. Recently, the effective potential 
for this condensate in the Curci-Ferrari gauge has been constructed in 5 by combining the 
algebraic renormalization 6 with the local composite operators technique [HI El) resulting in 
a dynamical mass generation. In this formalism, an essential step is the renormalizability of 
the local composite operator related to the condensate, which is fundamental to obtaining its 
anomalous dimension. It is worth mentioning that the anomalous dimension of the gluon-ghost 
operator O = {^A^A^ + ac a c a ^j in the Curci-Ferrari gauge, and thus of the gluon operator A 2 ^ 

in the Landau gauge, has been computed to three loops in the MS renormalization scheme, jHj. 
In addition, it has been proven by using BRST Ward identities that the anomalous dimension 
7^2 (a) of the operator A 2 ^ in the Landau gauge is not an independent parameter, being expressed 
as a combination of the gauge beta function, /3(a), and of the anomalous dimension, 7a(o), of 
the gauge field, according to the relation 

T*.(a) = " (^T + 7A(a)), " = ^ • (1-1) 

The aim of this paper is to extend the analysis of [H] to the Curci-Ferrari and maximal abelian 
gauges. We shall prove that the operator {^A aiJ, A^ l + aTfc 11 ^ is multiplicatively renormalizable 
to all orders of perturbation theory. Furthermore, as in the case of the Landau gauge, its 
anomalous dimension 70(a) is not an independent parameter of the theory, being given in fact 
by a general relationship valid in the Landau, Curci-Ferrari and maximal abelian gauges, which 
is 

7o(a) = - 2( 7c (a) + 7gc 2(a)) , (1.2) 

where 7 c (a), 7 9C 2(a) are the anomalous dimensions of the Faddeev- Popov ghost c a and of the 
composite operator ^gf abc c b c c corresponding to the BRST variation of c a . In other words 
sc a = ^gf abc c b c c where s is the BRST operator. 

The paper is organized as follows. In section 2, the renormalization of the dimension two 
operator ^A^A® + ac a c a ^J is considered in detail, by taking the Curci-Ferrari gauge as an 

In the case of the maximal abelian gauge, the color index a runs only over the N(N — 1) off-diagonal 
components. 
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example and the relationship (|1.2j) is derived. In section 3, we verify the relation we obtain 
between the anomalous dimensions explicitly at three loops in the Curci-Ferrari gauge in the 
MS scheme. Section 4 is devoted to the Landau gauge, showing that the expression (jl,2j) 
reduces to that of (jl.lj) . In section 5 we shall analyse the maximal abelian gauge, where the 
results established in iQ] for the case of SU{2) will be recovered. Finally, in section 6 we present 
our conclusions. 



2 The gluon-ghost operator in Yang-Mills theories in the Curci- 
Ferrari gauge 

2.1 The Curci-Ferrari action 

We begin by reviewing the quantization of pure SU(N) Yang-Mills in the Curci-Ferrari gauge. 
The pure Yang- Mills action is given by 

Sym = ~ \Jd 4 xF a ^F« u , (2.3) 

with F£ v = d^Al - d u A a ^ + gf abc A^Al. The so called Curci-Ferrari gauge, HUH2, is defined 
by the following gauge fixing term 

% = J d 4 x (b a d^ a + ^b a b a + c a d^ (D^c) a 

ry gf abc b a c b c c --g 2 f abc c a c b f cmn c m c n ) , (2.4) 



2 

with 

(ZV) a = d,c a + gf hc A\c c . (2.5) 

In order to analyse the renormalization of the operator {^A a>J 'A a ^ + ac a c a ^j , we have to introduce 

it in the action by means of a set of external sources. Following it turns out that three external 
sources J, ry M and r M are required, so that 

Sj = J d 4 x (j {^A a »Al + ac a c a ^j + \j 2 ~ ^A^ - t» (d^c a ) c a 

- 9 -^r bc A1c b c^ (2.6) 

where £ is a dimensionless parameter, accounting for the divergences present in the vacuum 
Green function (^^A a ^A^ + ac a c a ^J (^A a ^A a ^ + oic a c a ^j y, which are proportional to J 2 [3]- 
The action (Sym + 5gf + <Sj) is invariant under the BRST transformations, which read 



sAl = - {D,c) a , sc a = y -f abc c b c c , 



9 tabcjy 
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si? = b a , sb a = -Jc a , 

ST » = - rf , sr f = d^J , 

sJ = . (2.7) 

Also, due to the introduction of the external sources J, rf 1 and it follows that the operator 
s is not strictly nilpotent, namely 

s 2 $ = o , ($ = A a », c a , J, ^) , 
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s 2 c a = - Jc a , 

s 2 b a = - Jp abc c b c c , 

s 2 r" = - d»J . (2.8) 

Therefore, setting 

s 2 = 5, (2.9) 

we have 5 (Sym + Sgf + 5j) = 0. The operator 5 is related to a global SL(2,R) symmetry [T3] . 
which is known to be present in the Landau, Curci- Ferrari and maximal abelian gauges, |14j . 
Finally, in order to express the BRST and 5 invariances in a functional way, we introduce the 
external action 

S cxt = J <fx (n atl sAl + L a sc a ) (2.10) 
= J d 4 x (- n a » (D^cf + L a |/ a6c c 6 c c ) , 

where fi a/i and L a are external sources invariant under both BRST and 5 transformations, 
coupled to the nonlinear variations of the fields A® and c a . It is worth noting that the source 
L a couples to the composite operator |/ a6c c b c c , thus defining its renormalization properties. It 
is easy to check that the complete classical action 

E = Sym + Sgf + 5j + Sext (2-11) 

is invariant under BRST and 5 transformations 

sE = , SH = . (2.12) 

When translated into functional form, the BRST and the 5 invariances give rise to the following 
Ward identities for the complete action E, namely 

• the Slavnov- Taylor identity 



5(E) = , (2.13) 

with 

5 S = d A x [ — — — — — + — — + b a — - Jc a — + d^J- rf — 2.14 

v ; J \8n a i*8A« 5L a Sc a 5c a 5b a 6r)» 5t^ J v ' 

The 5 Ward identity 

W(E) =0, (2.15) 



with 



W(E) = / d'x (-Jc^ - - • (2.16) 

y ' J V $c a SL a Sb a 5t^J y ' 
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2.2 The invariant counterterm and the renormalization constants 



We are now ready to analyse the structure of the most general local counterterm compatible 
with the identities ([2.1.3ft and (|2.15j) . Let us begin by displaying the quantum numbers of all 
fields and sources 
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3 



(2.17) 



In order to characterize the most general invariant counterterm which can be freely added to all 
orders of perturbation theory, we perturb the classical action E by adding an arbitrary integrated 
local polynomial £ count in the fields and external sources of dimension bounded by four and with 
zero ghost number, and we require that the perturbed action (E + eE count ) satisfies the same 
Ward identities and constraints as E to first order in the perturbation parameter e, which are 



S(E + eE count 
W f E + eE count 



+ 0(e 2 ) , 
+ 0(e 2 ) , 



(2.18) 



This amounts to imposing the following conditions on E 



count 



£ S E C 



0. 



(2.19) 



with 



d 4 x 



5E 



ffi j 5T, 5 5E 5 

+ TT^Tir— + — — - + 



5 At SCl^ SQw 5 At 



6L a 5c a 5c a 5L a 



Jc a - 



5b a 



(2.20) 



and 



d A x -Jc a 



<5E C 



5c a 



^<5E<5E count ^<5E5E count ^5E count \ ^ 



5L a 



5b a 5L a 



5t^ 



(2.21) 



Proceeding as in 15], it turns out that the most general local invariant counterterm compatible 
with the Ward identities (|2.13j) and (|2.15j) contains six independent parameters denoted by a, 
0,2, as, 04 and CI5, and is given by 



^count 



(fx 



ZpWF^ + ( fl3 _ a A ){D^) a A a v + ^b a b a + a 2 b a d»Al 

+ (aa - a 3 )c a d 2 c a + (04 - a 2 )gf abc c a d» (c b A^) 

+ ( aaA ~ ai \ f^% a c b c c + (aa 4 - y ) f abc -^c b f cmn c m c n 
+ a^d^c a + a±gr bc ^A\c c - ^-L a gf abc c b c c 

+ ( ° 2 ± Q3) JA^A^ + ai Jc a c a + a 5 | J 2 - a 2 r/M£c a 

+ (02 - as^c^c 8 + (a 4 - a 2 )f r^/ afec A»c 6 c ( 



(2.22) 



It therefore remains to discuss the stability of the classical action. In other words to check that 
E co can be reabsorbed in the classical action E by means of a multiplicative renormalization 
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of the coupling constant g, the parameters a and £, the fields {<p = A, c, c, b} and the sources 
{$> = J, r, L, 0}, namely 

£(<?, e, a, *) + eS count = S( 50) £o, «o, <fo, *o) + 0(e 2 ) , (2.23) 
with the bare fields and parameters defined as 



a 

0/( 


= 


o a 

> s %t 


= Z n^1 


, 7"0/z = 




c 


= zyv 


ja 


= Z L L a 


, 9o = 


Z g9 , 


c 


= zJ'V 


3 <^0 


= ZjJ 


, a = 


Z a a , 


°0 


= Z 6 1/: V , 






, £o = 





(2.24) 

The parameters a, ai, a?,, a 3 , 04 and turn out to be related to the renormalization of the 
gauge coupling constant g, of the gauge parameter a, and of L a , c a , A^, and £ respectively, 
according to 

Z g = l + e a 



2 



Z a = 1 + £ I — — 0" — 2a2 + 2a 3 — 204 
Z L = 1 + e ( -a 2 - - + a 3 - a 4 



zy 2 = l + e 



-a 3 + 02 



2 



Z A 7 = l + e(-- + a 3 -a 4 



1/2 
4 

Z { = l + e(a 5 -2a-2a 2 + 2a 3 - 4a 4 ) . (2.25) 



Concerning the other fields and the sources f2°, 77^, and r^, it can be verified that they are 
renormalized as 

Zc = Z c , Z^ 2 = Z L X , Z^ = z^z^^z*/ 2 

Z„ = Z L 1 Z 6 7 1/2 , Z T = 1 . (2.26) 

Finally, for the source J, one has 

Zj = ZfZ-\ (2.27) 

from which it follows that 

70(a) = - 2( 7c (a) + 7^ (a)) , (2.28) 

where 7 c (a) and 7 9C 2(a) are the anomalous dimensions of the Faddeev-Popov ghost c a and of 
the composite operator ^f abc c b c c , defined as 



7 c (a) = /i<9 M lnZ c 1/2 l gc 2{a) = ^9 M lnZ L 70(a) = In Zj 

= (idplnZ- 1 7 a(a) = ^lnZ- 1 (2.29) 



a 

where ^ is the renormalization scale 



Therefore, we have provided a purely algebraic proof of the multiplicative renormalizability 
of the gluon-ghost operator to all orders of perturbation theory. In particular, we have been able 
to show, as explicitly exhibited in (|2.28[) . that the anomalous dimension of (j^A^A® + ac a c a ^j 
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is not an independent parameter of the theory, being given by a combination of the anomalous 
dimensions j c ( a ) and r y gc 2(a). It is worth mentioning that it has also been proven, ^3], for the 
Curci-Ferrari gauge that the anomalous dimension of the ghost operators ^f abc c b c c , ^f abc c b c c 
and ^f abc c b cf are the same. 

Although we did not consider matter fields in the previous analysis, it can be checked that 
the renormalizability of O and the relation (j2.28|) remain unchanged if matter fields are included. 



3 Three loop verification 



In this section, we will explicitly verify the relation (|2.28|) up to three loop order in the Curci- 
Ferrari gauge in a particular renormalization scheme, MS. The values for the /3-function and 
the anomalous dimensions of the gluon, ghost, the operator O and the gauge parameter a have 
already been calculated in the presence of matter fields in [Sj. For completeness we note that 
for an arbitrary colour group these are 



(3(a) 



7a(o) 



7c (a) 



70 (a) 



7a (a) 



^-C 2 A - AC F T F N f - ^-C A T F N f 



20 
3 



+ 



2 7\r2 



2830C%T F N f - 2857C1 + 1230C A C F T F N f - 2>IQC A T F N } 



W8C F T F N f -264C F T F Nj? — + 0(a 



54 



(3.30) 



[(3a - 13) C A + 8T F N f ] 



+ 
+ 



(a 2 + 11a - 59) C\ + A0C A T F N f + 32C F T F N f ] ^ 

54a 3 + 909a 2 + (6012 + 864C(3))a + 648C(3) - 39860) C\ 

- (2304a + 20736C(3) - 58304) C 2 A T F N f + (27648C(3) + 320) C A C F T F N f 



- 9728C A T F Nf - 2304C F T F N f - 5632C F T F Nf 
(a - 3)C A ~ + [(3a 2 - 3a - 95) C\ + 40C A T F N f 



1152 



+ 0(a 4 



(3.31) 



48 



+ 



(162a 3 + 1485a 2 + (3672 - 2592C(3))a - (1944((3) + 63268)) C\ 



(6048a - 62208C(3) - 6208) C 2 A T F N } - (82944C(3) - 77760) C A C F T F N f 
3 

+ 0(a 4 ) , (3.32) 



+ m\QC A T F Nf 



6912 



[16T F N f + (3a - 35)Ca] 



6 



a 

24 



280C A T F Nj + (3a 2 + 33a - U9)C\ + 192C F T F N f 

((2592a + 1944)C(3) + 162a 3 + 2727a 2 + 18036a - 302428) C\ 

- (62208C(3) + 6912a - 356032) C 2 A T F N f + (82944<C(3) + 79680) C A C F T F N f 

- mm8C A T F Nf - U82AC 2 F T F N f - 33792C F T F N^ 
2 

16 



3456 



+ 0(a 4 



(3.33) 



C A + (a + 5)C 2 J 



7 



+ 3Cl [(a 2 + 13a + 67a) C A - 40T F N f 



128 



+ 0(a 4 



(3.34) 



where the anomalous dimension of O in our conventions is given by (—4) times the result quoted 
in 0. The group Casimirs are tr (T a T 6 ) = T F 5 ab , T a T a = C F I, f acd f bcd = 5 ab C A , N f is the 
number of quark flavours and C( n ) is t ne Riemann zeta function. Our definition here of 7a (a), 
which denotes the running of a, differs from that of |Hj due to a different definition of Z a . For 
computational reasons, it turns out to be more convenient to consider the renormalization of 



the ghost operators f c c c , f c c c 



and / fl V? instead of gf abc c b c c , gf abc c b c c 



and gf abc c b c c 



respectively. We note that we have first verified that to three loops the anomalous dimension of 
each of the three operators are in fact equal, in agreement with ^Sj. Accordingly, we find 



7 c2 (a) 



-C A a + [(18a + 95) Cj - 40C A T F iV) 



a 

24 



+ 



(621a 2 + (7182 + 2592C(3))a + (1944C(3) + 63268)) C\ 
- (432a + 62208C(3) + 6208) C 2 A T F N f + (82944C(3) - 77760) C A C F T F N f 



- %%mc A T F Nf 



3456 



+ 0(a 4 



(3.35) 



We have deduced this result using the Mincer package, ^)], written in Form, ^7], where 
the Feynman diagrams are generated in Form input format by Qgraf, ^5j. For instance, for 
f abc c b c c there are 529 diagrams to determine at three loops and 376 for the operator f abc c b c c 
where each is inserted in the appropriate ghost two-point function. The same Form converter 
functions of |Hj were used here. Since the operator f abc (^c c has the same ghost structure as 
the operator c a c a , we were able merely to replace the Feynman rule for the operator insertion 
of c a c a in the ghost two-point function with the new operator and use the same routine which 
determined 70(a) in |H][. However, as f abc c b c c has a different structure we had to generate a 
new Qgraf set of diagrams to renormalize this operator. That the anomalous dimensions of 
both operators emerged as equivalent at three loops for all a provides a strong check on our 
programming as well as justifying the general result of section 2. Now, taking into account the 
extra factor g, the anomalous dimension ^ gc 2 (a) is found to be 



lqA a ) 



[8T F N f 
+ 



13CU] - + 

D 



(6a - 59) C\ + 40C A T F N f + 32C F T F N f — 
u -18 

(207a 2 + (2394 + 864C(3))a + (648C(3) - 39860)) C\ 
- (144a + 20736C(3) - 58304) C 2 A T F N f + (27648C(3) + 320) C A C F T F N f 



9728C A T F N? - 230AC F T F N f - 56326^7^ 



1152 



0(a 4 



(3.36) 



It is then easily checked from the expressions (|3.33|) . (|3.34|) and (|3.36|) that, up to three loop 
order, 

lo(a) = - 2( 7c (o) + l g Aa)) ■ (3.37) 

It is worth mentioning that the renormalizability of the operator O was already discussed in 
|19j from the viewpoint of the massive Curci-Ferrari model. Whilst the relation (|2.28|) was not 
explicitly given in ^H^, it is possible to obtain the relation from that analysis. Although the 
relation (|2.28f) has been established in the case of the Curci-Ferrari gauge, it expresses a general 
property of the gluon-ghost operator which remains valid also in the Landau and maximal 
abelian gauges, as will be shown in the following sections. 
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4 The Landau gauge 



The Landau gauge is a particular case of the Curci-Ferrari gauge, corresponding to a = 0. The 
Landau gauge is known to possess further additional Ward identities 01201) implying that the 
renormalization constants Zi, and Z c can be expressed in terms of Z g and Za, according to [0] 

Z L = Z\ 12 , Z c = Z~ l Z2 1/2 . (4.38) 

Therefore, it follows that (|2.27|) reduces to 

Zj = Z g Z A 1/2 (4.39) 

from which the expression is recovered, providing a nontrivial check of the validity of the 
general relationship (|1.2j) . 

As another internal check of our computations, we note that we should also find in the 
Landau gauge that r ) gc 2(a) = 7.4(0), as is obvious from ()4.38j) . It can indeed be checked from 
(EHTTT) and (|3~3H|) that 

V 2 ( a )| Q=0 = 74(a) | a=0 . (4.40) 



5 The maximal abelian gauge 



As is well known, the maximal abelian gauge is a nonlinear partial gauge fixing allowing for a 
residual U(l) N ~ l local invariance |1U1 1211 l2*2l I23j . In the following, a Landau type gauge fixing 
will be assumed for this local residual invariance. The Slavnov- Taylor and the S Ward identities 
H2.13j) and (|2.15|) can be straightforwardly generalized to this case. It is useful to recall that the 
gauge field is now decomposed into its off-diagonal and diagonal components 

A a T a = A irpi + A a T a ^ 

where the index i labels the N — 1 generators T 1 of the Cartan subalgebra of SU(N). The 
remaining N(N — 1) off-diagonal generators T a will be labelled by the index a. Accordingly, 
for the Faddeev-Popov ghost c a we have 

c a T a = c i T i + c a T a , (5.42) 

with 

sc a = gf^cPj + ^f a ^c c', (5.43) 
2 J 

Also, the group index of the gluon-ghost operator runs only over the off-diagonal components, 
namely 

Omag = (^A a »A« + ac a c a ^ . (5.44) 

Denoting respectively by Z and Z the renormalization factors of the off-diagonal and diagonal 
components of the fields, it follows that, according to the relationship (|2.28|) . the output of the 
Slavnov- Taylor and 5 Ward identities gives 

7ew( a ) = - 2 (7c Q (a) +7 9C 2(a)) , (5.45) 
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where 7 C "(a) and 7 9C 2(a) are the anomalous dimensions of the off-diagonal ghost c a and of 

the composite operator (^g f c /3 c 1 + |/ a ^ 7 c^c 7 ^ which corresponds to the BRST variation of 

c a . Moreover, as shown in [22], the use of the Landau gauge for the local residual U(1) N ~ 1 
invariance allows for a further Ward identity. This identity, called the diagonal ghost Ward 
identity in [22], implies that the anomalous dimension j gc 2(a) can be expressed as 

IgAa) = - 7c Q (o) - 7 C * (o) , (5.46) 

a a 

where J c i(a) is the anomalous dimension of the diagonal ghost c\ Therefore, for the expression 
of 7e> MAG (a) we obtain 

7o MAG (a) = - 2^ - 7ci (a)) , (5.47) 

a result which is in complete agreement with that already found in ^U] for the case of SU(2). 
Finally, it is worth mentioning that the anomalous dimensions of the diagonal and off-diagonal 
components of the fields have been computed at one-loop order in so that (|2.28j) gives 

explicit knowledge of the one-loop anomalous dimension of the gluon-ghost operator in the 
maximal abelian gauge. 



6 Conclusion 

We have shown that the mass dimension two gluon-ghost operator O = ^A^A^" 1 + ac a c a is 
multiplicatively renormalizable in the Landau, Curci-Ferrari and maximal abelian gauges. Fur- 
ther, we were able to establish a general relation between the anomalous dimension of O, the 
Faddeev-Popov ghost c a and the dimension two ghost operator gf abc c b c c , as expressed by the 
ea. (jl,2j) . This relation has been derived within the framework of the algebraic renormaliza- 
tion UJ], following from the Slavnov- Taylor identity (|2.13|) , As such, it extends to all orders of 
perturbation theory and is renormalization scheme independent, for any scheme preserving the 
Slavnov- Taylor identity. It has been explicitly verified up to three loops in the MS scheme in 
the Curci-Ferrari gauge. 

Furthermore, due to additional Ward identities that exist in the Landau gauge and in 
the MAG |22| . we were able to rewrite the relation (|1.2j) for the anomalous dimension for O 
in terms of the /3-function and the anomalous dimension of the gluon and/or ghost fields. In 
particular, concerning the maximal abelian gauge, it is worth underlining that the multiplicative 
renormalizability of the gluon-ghost operator, eq. (|5.47j) . is a necessary ingredient towards the 
construction of a renormalizable effective potential for studying the possible condensation of the 
gluon-ghost operator and the ensuing dynamical mass generation, as done in the Landau |Hj and 
Curci-Ferrari [3] gauges. 

As a final remark, we point out that, from the 3-loop expressions given in section 3, it is 
easily checked that the following relations holds in the Curci-Ferrari gauge: 

7o(o) = " + 7A(a)) , (6.48) 

7 9c2 (a) = 7x (o) - 2 7a (a) . (6.49) 

Up to now, we do not know if these relations are valid to all orders. They do not follow from 
the Slavnov- Taylor identity ()2.13|) . Nevertheless, although eas. (|6.48|) . (|6.49l) have been obtained 
in a particular renormalization scheme, i.e. the MS scheme, it could be interesting to search 
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for additional Ward identities in the Curci-Ferrari gauge which, as in the case of the Landau 
gauge 9 , could allow for a purely algebraic proof of eas. (|6.48|) . (|6.49l) . Notice in fact that, when 
a = 0, eq. ()6,48j) yields the anomalous dimension of the composite operator in the Landau 
gauge. Also, ea. (|6.49|) reduces to the relation (j4.4Uf) of the Landau gauge, since j a (a) = if a 
= 0. 
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